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ABSTRACT 

We calculate the angular correlation function of galaxies in the Two Micron All Sky Survey. We 
minimize the possible contamination by stars, dust, seeing and sky brightness by studying their 
cross correlation with galaxy density, and limiting the galaxy sample accordingly. We measure the 
correlation function at scales between V < < 18° using a half million galaxies. We find a best fit 
power law to the correlation function has a slope of — 0.79±0.02 and an amplitude at 1° of O.lOiO.Ol in 
the range 1' — 2.5°. However, there are statistically significant oscillations around this power law. The 
largest oscillation occurs at about 0.8 degrees, corresponding to 700/i~^ kpc at the median redshift of 
our survey, as expected in halo occupation distribution descriptions of galaxy clustering. In addition, 
there is a break in the power-law shape of the correlation function at 9 > 2.5°. Our results are in 
good agreement with other measurements of the angular correlation function. 

We invert the angular correlation function using Singular Value Decomposition to measure the three- 
dimensional power spectrum and find that it too is in good agreement with previous measurements. A 
dip seen in the power spectrum at small wavenumber k is statistically consistent with CDM-type power 
spectra. A fit of CDM-type power spectra in the linear regime (fc < 0.15/i Mpc~^) give constraints 
of flmh = 0.13 ± 0.07 and erg = 1-0 ± 0.09 for a spectral index of 1.0. This suggest a /Cs-band linear 
bias of 1.1 ± 0.2. These measurements are in good agreement with other measurements of the power 
spectrum on linear scales. On small scales the power-law shape of our power spectrum is shallower 
than that derived for the Sloan Digital Sky Survey. This may imply a biasing for these different 
galaxies that could be either waveband or luminosity dependent. The power spectrum derived here in 
combination with the results from other surveys can be used to constrain models of galaxy formation. 
Subject headings: galaxies: clusters: general — galaxies: statistics 



1. INTRODUCTION 

Correlation statistics are an important method for re- 
lating galaxies to the underlying mass distribution. The 
angular correlation function, w{6), measures the pro- 
jected clustering of galaxies by comparing the distribu- 
tion of galaxy pairs relative to that of a random dis- 
tribution. While a less direct probe than the three- 
dimensional correlation function ^(r), the angular cor- 
relation function can be a powerful approach owing to 
the larger sizes of two-dimensional surveys. The an- 
gular correlation function for bright galaxie s has bee n 
measured for the Lick Survey (Groth & Pee bles! [19771) . 
the Automated Plate Measuring (APM) Galaxy Sur- 
vey (|Maddox et al. 1990, 19^), the Edinburgh Durham 



Southern (EDS) Galaxy Catalog llCollins et al 
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and the Muenster Red Sky Survev ijBoschar 
Most recently t he Sloan Digital Sky Survey (SDSS; 
lYork et 'aD '2000') has performed a very detailed analy- 
sis of the angu lar correlat ion function ( Co nnoUv et al.l 
|2P02; Scranto n et alJl2002D from t heir Early Data Re- 
lease (EDR; S toughton et al.ll2002(l . 

Additionally, numerous papers have been written on 
techniques for inverting the angular correlation function 
to determine th e full t hree-dimens ional power spectrum, 
P(k) fUimbe ^ [T95I ILucvI pTi Jgaugh & Efstathiou 
1993[lpodelson k Ga ztaiiagal2000tlEfstathiou & Moodv 
200 it lEisenstein fc Zaldarriaga 2001; Padilla & Baugh 



2003) 



In this paper we calculate the angular correlation func- 
tion for galaxies selected from the Two Micron All Sky 



Survey (2MASS; ISkrutskie et"an I1997|) . We then in- 
vert w{9) using Singular Value Decompos ition, as sug- 
gested bv lEisenstein fc Zaldarriagal l)2001(l . to measure 
the three-dimensional power spectrum. A preliminary 
analysis of w{9) and P{k) for the Second Incremental 
Data Release of the 2MASS catalog was performed by 
^Ugood et al. (200l|). Here we use the complete and fi- 
nal 2MASS catalog, which provides a more than two-fold 
increase in t he number of sources owing to the nearly full 
sky (> 98%. lJarrett et a l.""2000bl coverage, and removes 
th e significant complicat ion of coverage masking suffered 
bv lAUgood et alJ \200)i} . In addition, the final 2MASS 
pipeline processing has improvements in the flux mea- 
surements that affect the catalog completeness. 

The angular correlation function of the 2MASS galax- 
ies is interesting for a number of reasons. First, with 
2MASS we can characterize the large-scale clustering 
of galaxies in the near-infrared (NIR) using the Kg 
(2.15/^m) passband. Galaxy Kg light is 5-10 times less 
sensitive to dust and stellar populations than _B-band 
light, providing a more uniform survey of the galaxy pop- 
ulation. Because the if^-band most closely measures stel- 
lar mass, it is possible that the Kg selected correlation 
function is a better probe of the dark matter power spec- 
trum than B and i?-band selected measurements. The 
variations of correlation functions in different bands con- 
tain information about how galaxy properties are related 
to the underlying dark matter distribution and, there- 
fore, to their formation and evolution. 

Secondly, 2MASS has full sky coverage extending out 



2 



to a median redshift of z = 0.07 and, hence, measures 
the power spectrum of our local Universe. Knowledge 
of our local "cosmography" , and how it may differ from 
other regions of the Universe, is relevant for many cosmo- 
logical tests. Thus the angular correlation function and 
power spectrum from 2MASS are valuable tools for un- 
derstanding galaxy formation and cosmology especially 
when used in comparison with measurements from other 
surveys. 

Our paper largely follows the treatment of th e 
Sloan EDR p erformed bv iScranton et all l)2002j) . 
iConnoUv et"aLri|2002f) and lDodelson et al.l l|200^ and is 

outlined as follows: In 33 we describe our galaxy sample 
selection from 2MASS. In []21we analyze the importance 
of systematic errors by studying their cross correlation 
with galaxy density. In 21 '^fi calculate the correlation 
function and discuss its errors, and in SjSlwe invert w{0) 
to measure the three-dimensional power spectrum, P{k). 
We conclude in [J^l 

2. 2MASS SELECTED GALAXIES 

To measure the angular correlation function accu- 
rately, it is critically important to fully understand the 
reliability and completeness of the galaxy sample. We 
select galaxies from the 2M ASS extended source catalog 
fXSC: IJarrett et al.ll2000b(l . which contains over 1.1 mil- 
hon extended objects brighter than Kg — 14 mag.^ The 
detection of galaxies is limited predominantly by confu- 
sion noise from stars, whose number density increases 
exponentially towards |6| = 0°. The XSC is mostly 
galaxies at |6| > 20° (> 98%), with an increasing stellar 
mixture (- 10%) at lower latitudes of 5° < |&| < 20°. 
Within the Galactic plane, i.e. |&| < 5°, there is an 
additional contamination by artifacts (10 — 20%) and a 
variety of Milky Way extended sources (~ 40%) includ- 
ing globular and open clusters, planetary nebulae, H II 
regions, you ng stellar objects, ne bulae, and giant molec- 
ular clouds ijJarrett et al.ll2000a|) . For our final analysis, 
we employ a latitude cut of |6| > 20° to remove the strong 
contamination of w{9) from stars at lower latitudes (see 
§3.1). Our cut reduces the Galactic contaminant sources 
to ^ 2% for our galaxy sample. 

The XSC also contains a fraction of spurious sources 
comprised of multiple star systems. Multiple star sys- 
tems become an increasingly important source of contam- 
ination as the stellar density, ristar rises, where rigtar is the 
number of Kg < 14 mag stars per sq. degree calculated 
in a co-add which is 8.5' x 16'. Specifically, the reliability 
of separating stars from extended sources drops rapidly 
from 95% at |&I = 10° (ristar ~ 5,000 deg"^), to < 65% 
at |6| < 5° (nstar > lO'^deg"^) l|Jarrett et al.ll2000a|) . To 
remove sources that are certainly unreliable, we use a cut 
at TT-star < 10^ dcg"^ during our initial galaxy selection 
from the XSC. The stellar density is this high over 5.5% 
of the sky so our galaxy catalog only covers 94.5% of the 
sky At high latitudes (|&| > 20°, ristar < l,250deg"^) 
the catalog is > 98% reliable for sources with Kg < 13.5 
ijJarrett et al.ll2000b(l . The XSC contains a smaU frac- 
tion of artifacts, e.g. diffraction spikes, meteor streaks, 
infrared airglow, at all latitudes. The XSC processing 

^ Completeness in the XSC is determined in terms of magnitudes 
measured inside the 20 mag arcsec"-^ eUiptical isophote (k_m_k20fe 
in the database). 
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Fig. 1. — The density of 2MASS galaxies on the sky in bins of 
Galactic longitude (/), and latitude (b). The horizontal error bars 
represent the width of each bin, while the Poisson uncertainties 
are given by vertical error bars. The dotted lines in the right-hand 
panel delineate the \b\ < 15° region corresponding to the Galactic 
plane. The dashed line shows the average galaxy density of 20.2 
deg~^. There is a significant variation between regions but there 
is no trend with longitude or latitude. 



removes and/or flags most artificial sources, as described 
fully in Jarrctt et al. (2000b). We use the XSC confusion 
flag (cc_flag) to remove sources identified as artifacts. We 
also remove a small number of bright {Kg < 12 mag) 
sources 2000) with dust corrected J — Kg colors less 
than 0.7 and greater than 1.4, which were determined to 
be no n-extragalactic extended sources (sec Mailer et aQ 
|2Q03). Finally, selecting in the Kg-hand minimizes the 
inclusion of false sources caused by infrared airglow. 

The XSC meets the original 2MASS science require- 
ments: greater than 90% completeness for extended 
sources (extragalactic and Galactic) with Kg < 13.5, and 
free fr om stellar confusion for |6| > 20° ijjarrett et alJ 
I2000b(l . In practice, the 2MASS completeness IS a mea- 
sure of the fraction of extended sources of a given magni- 
tude that are actually detected. Huchra & Mader (2000, 
private communication)^ show that the XSC is 99% com- 
plete for \b\ > 30° sources in the range 12.0 < Kg < 13.7 
mag, using log A^-log S and V/Vmax tests. At magnitudes 
brighter than Kg — 11.5, they deter mine that greater 
than 95% of the known gal axies in iZwickv fc Kowall 
ll96g) are found in the XSC. I.Tarrett et all l|2000aD find 
the XSC completeness remains at 95% for Kg < 13.0 
sources well into the Ga lactic plane (|6| > 5°). Addi- 
tionally, iBell et all 1I200I cross correlate the XSC with 
the complete (> 99% lStrauss et alJl2002D . spectroscopic 

flaxy sample drawn from the SDSS EDR. Bell et aQ 
0037 find that in the large (414 square degree), high 
latitude (|6| > 30°) region of the EDR, the XSC misses 
only 2.5% of the known galaxy population down to 
Kg < 13.57 in dust correc t ed Kr on magnitudes. 

2MASS employs E roni l|1980l) magnitudes to obtain 
the total flux of each galaxy. The 2MASS-defined Kron 
magnitudes (k_m_fe in the XSC) are elliptical aperture 
magnitudes with semi-major axes equal to 2.5 times the 

^ |cf a-www ■ harvard ■ 6du/~huchra/2mass/verif y ■ htm| 
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first moment of tlie briglitness distributions for each 
source. Tlie first moment calculation, the Kron ra- 
dius, is computed to a radius that is five times the 20 
mag arcsec"^ isophote. 2MASS limits the Kron radius 
to a 5" minimum (owing to the PSF; see the 2MASS 
Explanatory Supplement for details). Because of the 
short exposure times (7.8 seconds with a 1.3-m telescope, 
iSkrutskie et al.lll997|) of 2MASS, the Kron magnitudes 
underestimate the true tota l flux systemati cally by ap- 
proximately 0.1 magnitudes fBell et al.*2003f). At a mag- 
nitude limit slightly fainter than the 2MASS science re- 
quirement. Kg ^ 13.57, the XSC galaxy completeness is 
quite good, primarily missing'^ only low surface bright- 
ness (1.5%) and distant (1.0%) gal axies both of whi ch are 
near the Kron = 13.57 limit l)Bell et al.l 120031. We 
adopt Kron magnitudes for this paper and henceforth all 
cited magnitudes will be Kron magnitudes. The small 
fraction of low surface brightness galaxies miss e d in th e 
XSC confirms the expectation of I.Tarrett et al.l l)2000hD . 
Moreover, the automated processing of the XSC pro- 
duces systematically incomplete photometry for galaxies 
larger t han 50" due to the typical 2MASS scan width 
of 8.5' l|.Ta,rrett et a,].ll2000hf> . We include the 2MASS 
Large Galaxy Atlas sample (540 galaxies identified in 
the XSC database with cc_flag='Z') from ijJarrett et al.l 
I2OO3') ■ which has been assembled to account for most of 
the scan size photometric incompleteness. 

We select initially all Kg < 14 mag objects from 
the XSC with ngtar < lO'^deg"^ . We apply Galactic 
foreground dust corrections from 'Schle gel et al.l l|1998f) 
to each galaxy in the Kg < 14 catalog and cut at 
Ks = 13.57 mag, resulting in 775, 562 galaxies. The 
dust correction increases the number of galaxies brighter 
than Ks = 13.57 by 10%. We plot the galaxy density as a 
function of longitude and latitude in Figure ^ illustrat- 
ing that the dust-corrected number densities are fairly 
uniform for |&| > 15. Ultimately, our final sample se- 
lection includes 501,578 galaxies with Kg < 13.57 mag, 
\b\ > 20°, and Galactic dust corrections less than 0.05 
magnitudes in i^^-band (as described in §3.2). We sum- 
marize these sample cuts leading to our final selection in 
Tabled 

The redshifts of many 2MASS galaxies have been mea- 
sured as part of other redshift surveys. Large red- 
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1,295,895 
775,562 
688,464 
618,333 
548,353 
480,794 
501,578 
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< 13.57 

< 13.57 

< 13.57 

< 13.57 

< 13.57 

< 13. 57=^ 



> 

> 

> 10 

> 15 

> 20 

> 25 

> 20 



''Only sources with extinction correction < 0.05 mag. 

TABLE 1 

Galaxy Sample Selection. The number of galaxies in our 
2MASS catalog is given as a function of Ks magnitude cut 
AND Galactic latitude cut. See ^ 



^ Some of these missed distant objects can be found in the 
2MASS point source catalog. 
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Fig. 2. — Galaxy-star cross-correlations (open diamonds) for 
cuts in Galactic latitude of |b| > 10°, 15°, 20° and 25°. Error bars 
displayed for the 25° cut are representative of the other cuts. The 
auto-correlation is shown as naked error bars for comparison. We 
require a latitude cut of |6| > 20° to be free of stellar contaminants 
out to separation angles of 9 < 15°. 



shift samples provide the redshift distribution dN/dz of 
2MASS galaxies, and allow calculations of the NIR lumi- 
nosity function jKochanek et al.ll20?nt iCole et al.lf20nit 
iBell et ;ini2003|) . For computing the three-dimensional 
power spectrum (see Et, we will use redshifts from 
the SDSS EDR ( S toughton eTal] 12002') and the 2dF 
Galaxy Redshift Su rvey (2dFGRS) 100k Data Release 
llColless et alJl20?)l . 

3. CROSS-CORRELATION FUNCTIONS 



FoUowing lScranton et al. (2002, hereafter Sc02) we an- 
alyze the cross correlation between the galaxy number 
counts and possible sources of systematic error. This 
is an effective method for identifying errors, quantifying 
their magnitude, and determining the selection limits re- 
quired to reduce them to manageable levels. Below we 
analyze the cross-correlation signal for four possible con- 
taminants: stars, dust, seeing, and sky brightness. 

To perform the cross correlation, we pixelize the galaxy 
sample on a regular grid to determine local galaxy 
densities. We u se the HEALPix ^ software package 
ijGrski et al.lll999() to create equal area pixels on the sky. 
We use a Nside parameter of 64 generating Nc = 49, 152 
cells each with an area of 0.84deg^. This yields an av- 
erage of 15.6 galaxies per cell so very few cells contain 
no galaxies. We calculate the average value of each con- 
taminant for a cell using the individual measurements for 
each galaxy within the cell. When making a latitude cut, 
we include all cells whose centers are above the given lat- 
itude including galaxies that might be below the cut. To 
calculate errors, we subsample the data, dividing the sky 
into 48 tiles, with each tile containing 1024 cells, making 
them roughly 29.3° on a side. Because we will perform 
a latitude cut and 8 of the tiles lie along the galactic 
plane we will end up only using 40 tiles for our analy- 
sis. For each tile we calculate the fractional galaxy and 

* http:/ /www. eso.org/science/healpix/ 
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Fig. 3. — Galaxy-extinction cross-correlations (open triangles) 
for three cuts in the maximum allowed value for the extinction 
given by A(Ks) = 0.367 X E{B - V), with E{B - V) from 
ISchleEfel et al.l 119981) . Error bars are shown only for the lowest 
cut and are representative. The auto-correlation is shown as naked 
error bars for comparison. A cut of A{Ks) < 0.05 mag reduces 
the cross-correlation signal to less then the auto-correlation for 
e < 10°. 



contaminant over density in a cell i by 

n ^ — m 9 



n9 
xf — x'^ 



(1) 



where the averages are calculated for each tile. The cross- 
correlation function, Wgc{0), is then 







1-3 



(2) 



where Qfj is unity if the separation between cells i and 
j is within angular bin 6^ and zero otherwise. The error 
on the mean can be computed by 

1 ^* 

[^^gMf = J^H ^^^sciO) - Wg,A0)f , (3) 
* i=l 

where Nt = 40 is the total number of tiles. We use 
this cross-correlation function to quantify the systematic 
contributions of contaminants in the subsections below. 

3.1. Stars 

During the 2MASS pipeline processing, the stellar den- 
sity ristar — the number of Kg < 14 mag stars per square 
degree — is determined for each extended source on a co- 
add which is 8.5' x 16' in size. In regions of high stellar 
density, stars may be mistakenly identified as extended 
sources. For 95% of the sky, ristar provides an excellent 
measure of confusion noise in the XSC. The stellar den- 
sity saturates in 2MA SS at values of ristar > 3 x 10** deg^^ 
l|.Tarrett et al.ll200flal) . occurring near the Galactic center 
(l&l < 5°, 350° < I < 10°). RecaU that for our galaxy 
sample we make a more conservative stellar density cut 
of nstar < lO^'deg"^. 

To examine the importance of this contaminant wc 
look at the cross-correlation of the galaxy density and 



Fig. 4. — The galaxy-seeing cross-correlation (open squares) 
as a function of angle. The auto-correlation can not be plotted 
because it is > 30 times larger. Seeing variations are not a source 
of systematic error. 



the stellar density. Because ristar is a strong function 
of Galactic latitude we measure the cross-correlation for 
different Galactic latitude cuts and present the results 
in Figure [3 Including XSC objects with low latitudes 
(l&l < 15°) results in a higher amplitude galaxy-star 
cross-correlation at all angular scales than the galaxy 
auto-correlation. That this cross-correlation has a differ- 
ent slope than the galaxy auto-correlation suggests that 
the signal is coming from the stellar auto-correlation as 
multiple star systems are mistakenly identified as galax- 
ies. However, if the latitude cut is increased to |5| > 20°, 
then the auto-correlation signal is not contaminated sig- 
nificantly by stars at angular separations oi 9 < 15°. 
Thus, we adopt a |6| > 20° cut throughout this paper. 
Making this cut reduces the area of the survey to 67% 
of the sky or 28,960 square degrees. 548,353 galaxies 
survive this latitude cut. 

3.2. Dust 

Though we have attempted to correct our galaxy 
sample for th e effec ts of foreground extinction using 
iSchlegel et all l|1998D . dust may remain an important 
contaminant. There is a strong correlation between 
stellar density and dust extinction, both of which in- 
crease substantially towards the Galactic plane. Here 
we look at the cross-correlation between dust extinc- 
tion and galaxy density only at |6| > 20°, owing to the 
stellar contamination at lower latitudes. In Figure 13 
we show the galaxy-extinction cross-correlation for three 
cuts in magnitudes of Kg-hand dust extinction. To bring 
the galaxy-extinction cross-correlation below the auto- 
correlation we mask out those cells with i^g-band ex- 
tinction of A{Ks) = 0.367 X E{B - V) > 0.05, which 
amounts to 2,602 cells. This leaves and area of 23,345 
sq. degrees and 501, 578 galaxies. 

Our use of a dust correction to galaxy magnitudes 
does increase the dust cross-correlation. For example, 
the amplitude of the galaxy-extinction cross-correlation 
is reduced by about a third for a catalog without a dust 
correction. However, the dust correction increases the 
number of galaxies by 10%, and therefore the error in- 
troduced by the correction is a minor contribution to the 
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Fig. 5. — The galaxy-background (Xs) cross-correlation function 
as a function of angle. The background error correlates with galaxy 
density as expected, yet it is always significantly less than the auto- 
correlation (naked error bars). Sky background fluctuations are a 
negligible source of systematic error. 



systematic error budget. 

3.3. Seeing 

Great care has gone into ensuring that the final 2MASS 
release has a highly uniform seeing over the entire sur- 
vey area. The pipeline processing allowed tracking of the 
point-spread function (PSF) as a function of time.^ The 
mean radial size of the 2MASS PSF over a scan could 
vary over periods as short as a few minutes. Depend- 
ing on the local density of stars, the seeing was tracked 
on time-scales ranging from 2 — 30 seconds. Scans with 
poor seeing were reobserved under better conditions to 
maintain the seeing uniformity. For each source in the 
XSC, the pipeline measured the average radial extent of 
the PSF in each passband (called "ridge shape" and de- 
noted k_shO in the XSC for the if^-band), analogous to 
a full width at half maximum (FWHM) measurement 
of the varying PSF. The full details of the PSF charac- 
terizat ion and seeing tracking are given in l|Jarrett et al.l 
I200nbt) . 

We use k_shO to test for systematic errors in the an- 
gular correlation of galaxies caused by spatial variations 
in seeing. One might expect that observations made un- 
der poor seeing conditions would lead to more spurious 
galaxy detections, causing a correlation between seeing 
and galaxy density. For our sample the mean seeing is 
k_shO — 0.992 with 0.049 rms scatter, corresponding to 
a typical survey seeing of 2.1" FWHM with 5% uncer- 
tainty. Given this uniformity, we do not expect the see- 
ing to be a major source of systematic errors. One might 
want to check for seeing variations over the scan width of 
8.5'. However, even in the densest regions there are only 
a few 2MASS galaxies in such a small area so a cross- 
correlation cannot be performed. In Figure ^ we plot 
the galaxy-seeing cross-correlation, whose amplitude is 
more than 30 times less than the auto-correlation ampli- 
tude. Hence, seeing variations are a negligible source of 

^ 2MASS mapped the sky using overlapping scans roughly 8f5 
wide by 6° long. The scan direction followed the declination axis, 
with each scan covering 6° in approximately 6 minutes. 



systematic error. 

3.4. Sky Background 

In Ks-hand, the dominant source of background flux is 
thermal continuum emission from the atmosphere. While 
the Kg background is less severe and variable than the 
airgiow induced J and H-band backgrounds, the sky 
brightness at Kg is still quite high at 13.3 ± 0.3 mag 
arcsec"^. This is orders of magnitude brighter than the 
typical outer isophotes of 2MASS galaxies, which have 
UK > 19 mag arcsec"^. Furthermore, the Kg background 
can be variable and may still produce high-frequency fea- 
tures extending to tens of arcseconds. To mitigate these 
features, the 2MASS processing includes a sophisticated 
backgr ound fitting and remo val scheme that is described 
fully in lJarrett et aD l)2000bj) . 

The rms background error caused by Poisson noise, 
which depends on the variable sky brightness, might also 
produce sources of systematic error in calculating the an- 
gular correlation of 2MASS galaxies. The median back- 
ground local to each source is determined in each pass- 
band and given in the XSC database (k_back for Kg). For 
each source in our sample, we estimate the background 
error from co-added and resampled images 
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The factor of 6 represents the number of co-added frames 
comprising each source image. Each 2" pixel has been 
resampled {n^ = 4) into 1" pixels and smoothed with a 
kernel giving rise to the factor of fcf = 0.5853. In i^'s-band 
2MASS has an average gain of G = 8.2 electrons/DN and 
an average read noise of ctrn = 52 electrons (T. Jarrett, 
2001, private communication) ^ . 

A larger sky uncertainty increases the difficulty of 
galaxy detection and leads to a correlation between back- 
ground error and galaxy density. We find an average 
k_back of 475 DN with scatter 125 DN, corresponding to 
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Fig. 6. — The cross-correlations with all four investigated con- 
taminants using our cuts in Galactic latitude of fe > 20° and in 
extinction of Kg < 0.05 compared to the galaxy auto-correlation 
signal. On angular scales < 10° the auto-correlation amplitude 
is significantly higher then the cross-correlation with any of the 
contaminants. 
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Fig. 7. — The integral constraint correction, Aw{9), compared to 
the angular correlation function w{9) for the pixelized estimator. 
The integral constraint is small because it is suppressed by a factor 
where in this case Nc = 32, 512. Thus the integral constraint 
is not relevant for large surveys like 2MASS and SDSS except at 
the largest angular scales. 



a background uncertainty of 20.02 ± 0.13 mag arcsec ^. 
This a grees with the typical la error of 20.0 mag arcsec"'^ 
in Ks (|,Tarrett et ani2000hL l200l . The estimated back- 
ground error is quite uniform and the galaxy density and 
background error correlation is weak, as we show in Fig- 
ure The cross-correlation amplitude is always smaller 
than that of the galaxy auto-correlation and hence back- 
ground fluctuations are a negligible source of error in 
determining w{0). 

3.5. Summary of Systematic Effects 

Figure El shows the cross-correlations of the four pos- 
sible contaminants in comparison with With the 
Galactic latitude cut of |6| > 20° and the dust extinction 
cut of Kg < 0.05, we see that the contaminants are sig- 
nificantly below the auto-correlation signal for angular 
scales < 10°. We measure w{9) on larger scales, but 
caution that for 6 > 10° systematic errors might be as 
large as the signal. 

4. THE ANGULAR CORRELATION FUNCTION 

4.1. Estimators, Biases, and Errors 

Following Sc02 wc use two estimators for the angu- 
lar correlation function. The first, which we will refer 
to as the pair-estimator, compares positions of the ob- 
served galaxies (data) to the positions of random points 
by counting the number of pairs in an angular bin 9a nor- 
malized by the total number of pairs. We use one million 
randomly placed points for the comparison, avoiding the 
regions that were masked as described in §3. The an- 
gular correlation can then be comput ed using the the 
now standard estimator introduced bv iLandv fc Szaliivl 



w{0) 



(DD) - 2{DR) + (RR) 
(RR) 



(5) 



where DD, DR, and RR are the normalized number of 
data-data, data-random and random-random pairs in an 
angular bin 6. At large angular scales counting pairs be- 
comes computationally expensive, so instead we calculate 



Fig. 8. — The angular correlation function of 2MASS galax- 
ies. Bins evaluated using the pixelized estimator are shown with 
an X, otherwise the pairs estimator is used (crosses). The errors 
are based on Jackknife resampling. The best-fit power laws of 
the APM (long dash ; [Maddox et al. 1990) and SDSS (dot-dash; 
IConnoUv et ail2002H surveys are shown for comparison. The solid 



line is a power law with a slope of - 
at 1° of 0.10 ±0.01. 



-0.79 ± 0.02 and an amplitude 



the galaxy density in a cell as described in ^ (eq. |2Jl, 
which we refer to as the cell-estimator. Then the value 
of the angular correlation function is given by 



w{9) = 



(6) 



where again 0^ is unity if the angular separation of cell 
i and cell j is within 9a and zero otherwise. We use 
the pair-estimator for 9 < 5° and the cell-estimator for 
9 > 2° which gives three angular bins where we use both 
methods to check that they generate the same results. 
In the overlapping regions we will use the pair-estimator 
which should be more accurate because the data has not 
been smoothed. 

All estimates of w{9) are subject to a statistical bias 
referred to as the "integral constraint " ijPeeblesI Il980t 
IBernsteinll99llHui fc Gaztailagall999j) . This bias arises 
because the estimate of the mean number density in a 
given cell enters into the estimator of the angular corre- 
lation function nonlinearly. The integral constraint cor- 
rection. Aw (9), for the cell-estimator can be calculated 
to be (see Sc02), 



Aw{9p) = [1 + (2ci2 - 3)w{90)] ^ 



(7) 



where ci2 ~ 2 and the sum is over all cells. Since this 
value is suppressed by a factor of this bias is ex- 
tremely small for large surveys like 2MASS and SDSS. 
Figure [7| shows Aw{9) in comparison with 1^(9). One 
clearly sees that the integral constraint can safely be ig- 
nored at all angular scales where we can measure w{9). 

To calculate errors we use the Jackknife resampling 
method, the method that Sc02 found to be the best data- 
only method for estimating errors. We remove one tile 
of 1024 cells and then use the two estimators described 
above to determine w{9) on the remaining data set. We 
therefore make Nt — 40 measurements of the angular 
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Fig. 9. — The correlation matrix of the angular correlation func- 
tion of 2MASS galaxies. The correlation matrix is the covariance 
matrix normalized by its diagonal elements. One sees that for bins 
less than 0.06° the data is fairly independent. But at scales larger 
than this there is a strong correlation between different data bins. 



correlation function in 25 angular bins. To calculate er- 
rors we determine a covariance matrix out of the 40 mea- 
surements. The a/3 element of the covariance matrix is 
computed by 

N ^* 

* fc=0 

where W}^{9a) refers to the fcth measurement of the an- 
gular correlation function on the angular scale 9a- 

Because of the strong covariance between different an- 
gular bins, it is necessary to include it when fitting a 
model to the data. Thus we minimize the complete form 
of the statistic 

x^ = Y. [^(^") - ^™(^")] - w„,(e0)] (9) 

where W rn.(O) is a inodel of the angular correlation func- 
tion fsee iBernsteinllTOQ^ for a discussion). 

4.2. Results 

The angular correlation function of 2MASS galaxies 
is shown in Figure |H1 We fit a power law of form 
'w{9) — A9^^^ out to angular scales of 2.5° to the es- 
timated w{9). Note that the uncertainties of the power- 
law fit are substantially underestimated if one ignores the 
covariance between angular bins. Hence to fit a power- 
law form to our estimated angular correlation function 
we use Singular Value Decomposition (see i)5.2|l to re- 
move the oscillatory modes from the covariance matrix. 
We then find a best fit power law with an amplitude 
at 1° oi A = 0.10 ± 0.01 and 1 - 7 = -0.79 ± 0.02. 
However, we find that a power law is not a good fit to 
the angular correlation function (x'^/d.o.f. ~ 2.4 using 



Fig. 10. — The angular correlation function of 2MASS galax- 
ies multiplied by 9^ "^^. Small oscillations around the power-law 
form are apparent and are statistically significant. The horizontal 
dashed line with amplitude A = 0.10 is the best fit power law. 
Bins evaluated using the pixelized estimator are shown with an X, 
otherwise the pairs estimator is used (crosses). 



eq. El using the full covariance matrix (see discussion 
below). As shown in Table 12 the slope of w{9) from 
2MASS is slightly higher but in overall agree ment with 
the slopes foun d in AP M ("Madd ox et al.ll990() and SDSS 
EDR (Connoll ~t all 12002). The power-law fit to the 
2MASS data extends past 6* = 1 - 2°, where the APM 
and EDR angular correlation functions start to deviate 
from a power law. This is expected since the same phys- 
ical scale will correspond to a larger angular scale in the 
shallower 2MASS catalog. 

The amplitude of the 2MASS angular correlation func- 
tion is 3.6 (7.7) times larger than the amplitude found 
in the APM (SDSS EDR) surveys. We expected there 
would be an amplitude difference both because 2MASS 
galaxies are brighter than either the APM or SDSS EDR 
survey and because 2MASS galaxies have a shallower red- 
shift distribution. Assuming that the galaxies in all three 
surveys have the same 3-dimensional power spectrum, we 
can use Limber's equation (eq. Ill() to calculate that the 
change in amplitude in w{9) when going from the me- 
dian redshift of SDSS (0.17) to the median redshift of 
2MASS (0.074) is approximately a factor of 5. From Ta- 
ble 121 we see that the ratio of amplitudes is a factor of 
7.7, hence even after including the median redshift differ- 
ence the 2MASS w{9) amplitude is still 50% higher than 
SDSS. We attribute the remaining amplitude offset to the 
difference in luminosity of the galaxies being sampled. 
The average (r — Ks) color of a galaxy is approximately 
3.0 so our magnitude cut of Kg < 13.57 corresponds to 
r < 16.5. Therefore, the faintest 2MASS galaxies are 
two magnitudes brighter than the mean of the bright- 
est magnitude bin used in the SDSS study. More lumi- 
nous galaxies are more clustered than le ss luminous ones 
^Norberg et al.ll200Tt IZehavi et a,1jl2nn^. 

In Figure EH we plot the angular correlation function 
of 2MASS galaxies multiplied by 9°-'^^ to examine the 
deviations from a pure power law. The poor value of 
X^/d.o.f. for the power-law fit, the straight- horizontal 
line, may seem surprising since the magnitude of the 
deviations are similar to the error bars. However, the 
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Survey A I-7 Zm sample selection 

APM 0.028 ±0.003 -0.70 ±0.03 0.11 bj < 20 

SDSS 0.013 ±0.060 -0.74 ±0.04 0.17 18 < r, < 19 

2MASS 0.10 ±0.01 -0.79 ±0.02 0.074 < 13.5 

2MASS 0.08 ± 0.01 -0.79 ± 0.02 0.080 13.5 < Ks < 12.5 

2MASS 0.22 ± 0.02 -0.79 ± 0.02 0.047 12.5 < Ks < 11.5 

2MASS 0.48 ± 0.03 -0.48 ± 0.02 0.029 11.5 < < 10.5 

2MASS 0.15 ±0.01 -0.79 ±0.02 0.071 < 17.75 

2MASS 0.08 ± 0.006 -0.76 ± 0.02 0.079 fiR > 17.75 

TABLE 2 

Th e resul ts of fitting the the a ngular correlation function with a functional form A9^ ' FOR APM IMADDOX ET AL.I 
119961) . SDSS ICONNOLLY ET ALll2002fl and this work. We also SHOW THE FITTED PARAMETERS WHEN THE 2MASS SAMPLE IS DIVIDED 

BY EFFECTIVE SURFACE BRIGHTNESS AND INTO MAGNITUDE BINS. 



depicted error bars follow from the diagonal elements 
of the correlation matrix and do not describe indepen- 
dent variance. Owing to significant off-diagonal terms 
in the covariance matrix, the displacement of adjacent 
w{6) bins towards the power law fit gives a very large 
contribution to x^. The correlation matrix of w{9) is 
shown in Figure |2| The importance of off diagonal ele- 
ments is clearly seen in the figure. Thus, the oscillatory 
behavior scene in Figure 1101 is significant. These may 
be the result of wiggles in the power spectrum caused 
by baryon oscillations as recently detected in the 2dF- 
GRS power spectrum ijPercival et alJl200H) . The oscil- 
lations seen here, however, are at much smaller scales 
then those seen in the 2dFGRS and are in the nonlinear 
regime. Such oscillations are expected in halo occupa- 
tion distribution descriptio ns of galaxy cluster ing (Scliak 
l200(llBerlind fc Weinbergll2002t iZehavi et al.ll2004,1 and 
the dip centered around 0.8 degrees, corresponding to 
700/i~^ kpc at the median redshift of our survey, occurs 
at about the expected physical scale. 

We also calcrdate the angular correlation of 2MASS 
galaxies split into bins of AKg = 1 centered on ^^3=11, 
12, 13 and plot each w{6) in Figure 1771 The error bars 
at a given angular separation increase for brighter mag- 



nitudes owing to their lower numbers. Limber's equation 

1) Limbed|1953|) implies that apparently brighter galaxies 
should have parallel w(9) with higher amplitudes owing 
to their closer redshift distribution. Using the median 
redshift of the galaxies in each magnitude bin (see Ta- 
ble |2Jl we estimate that the amplitude ratio should be 
6 : 2.5 : 1, going from the brightest to the faintest bin. 
In Table |21 we present the parameters of the power-law 
fits to w{0) for each magnitude bin and confirm that the 
differences in amplitude can be fully explained by differ- 
ences in median redshift, consistent with true clustering. 

Because 2MASS is an all sky survey with 1% photo- 
metri c uniformity on large angular scales ijNikolaev et alJ 
|2000|) . we have two independent volumes, one for north- 
ern b > 20° and one for southern b < —20° Galactic 
latitudes. We can compute w{6) independently for the 
two hemispheres as a check on our procedure and to ex- 
amine the effects of cosmic variance, plotted in FigurcIT^ 
At angles less than 0.5° there is good agreement. How- 
ever, at larger angles the northern Galactic hemisphere 
shows more clustering, caused by real differences in the 
observed large scale structure b etween the northern and 
southern Galactic hemispheres ijMaller et al. 2003, Fig. 

2) . For small angular scales, however, this cosmic vari- 
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Fig. 11. — The angular correlation functions of 2MASS galax- 
ies divided into bins of AKs = 1 centered on i4's=ll,12,13. The 
error bars, which are smaller then the symbols, decrease with Ks 
magnitude since there are fewer bright galaxies. The slopes re- 
main the same while the amplitudes increase for brighter galaxies 
as expected. 



Fig. 12. — The angular correlation function of 2MASS galaxies 
for northern Galactic latitudes {b > 20° diamonds) and southern 
Galactic latitudes {b < —20° squares). For angular scales less then 
0.5° the two hemisphere yield consistent results. However, at larger 
angular scales there is stronger clustering in the northern Galactic 
hemisphere. 
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Fig. 13. — SDSS r-band concentration = r"9o/^50 versus effec- 
tive surface brightness fix = Kg -1-2. Slog j^q (27rr^). The horizontal 
dotted line indicates Cr = 2.6, chosen by SDSS to divide early 
from late-type galaxies. The vertical dashed line at fj,K < 17.75 
mag arcscc^^ shows our adopted morphology divider. 



ance is unimportant. 

One would like to investigate the dependence of w{9) 
on galaxy morphology. In the optical, galaxy morphology 
can be estimated using either galaxy color or concentra- 
tion. In the infrared, however, the {J— Kg) colors contain 
little information about the galaxy type, since all galax- 
ies have (J — Ks) colors that are t ightly peaked around 
1.0. Furthermore . IBell et al.l l)2003fl demonstrate that the 
2MASS concentration measurement is not very accurate 
for faint {Kg > 12) galaxies. Instead, as a measure of 
galaxy morphology we use effective surface brightness 



(10) 



where is the half-light radius of the galaxy, and Kg 
is the extinction-corrected Kron magnitude. For 6, 238 
2MASS galaxies also in the SDSS EDR we show in Figure 
^|that /iif is correlated with the SDSS optical concentra- 
tion Cr = rgo/r^o, where rgo and r^o are the radii within 
which 90% and 50% of the galaxy flux are contained, re- 
spectively. The SDSS collaboration has adopted Cr = 2.6 
to separate between early and late types in a rudimen- 
tary fashion (e.g. Stratcva et al. 2001; Kauffmann ct aj 
iOQJ). Using Cr > 2.6, we find that 2630/3203 (82%) of 
the high surface brightness {fix < 17.75 mag arcsec"^) 
galaxies are early-type. 

Dividing the sample at the median fj,K = 17.75, we 
plot the angular correlation functions of the two popula- 
tions in Figure^! The higher surface brightness galaxies 
are more clustered with an amplitude at 1° almost twice 
that of the lower surface brightness galaxies (see Table 
121). Unlike the difference in amplitude for the different 
magnitude bins, which are caused by differences in the 
median redshifts, our high and low surface brightness 
galaxies have nearly identical median redshifts of 0.079 
and 0.071, respectively. Hence, the higher fiK (early- 
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Fig. 14. — The angular correlation function of high and low 
surface brightness 2MASS galaxies. Galaxies with fiK < 17.75 
(diamonds) are more clustered and have a slightly steeper slope 
than galaxies with /ifc > 17.75 (triangles). 



type) galaxies truly are more clustered in three dimen- 
sions than the low (late-type) galaxies in our sam- 
ple. This increased clustering amplitude for early type 
galaxies has also bee n seen in the EDR t hree-dimensional 
correlation function l|Zehavi e t al.i200^ and is related to 
the morphology-density relation (Dressier 198(j). 

The slopes of the two angular correlation functions are 
marginally different; they disagree at the 90% confidence 
level. Such a disagreement could result from nonlocal bi- 
asin g or nonlocal causes for the morph ology-density rela,- 
tion ijNaravanan et alJ200r)HScherrer fc Weinberg 1998). 
It would be interesting to see if any galaxy formation sce- 
nario matches this slight disagreement. 

5. THE 3D POWER SPECTRUM 

In this section we invert the estimated angular corre- 
lation function to measure the three-dimensional power 
spectrum. The relationship between w{6) and P{k) can 
be expressed as 



w{e) 



dkkP{k)g{ke) 



(11) 



where the kernel g{k9) is given by 



gikd) 



1- r dzMk9Xiz))i^r-^F{z) (12) 



l)Limbedll95l iBaugh fc Efstathioulll993(l . In the above 
equation X{z) is the comoving distance to a redshift z, 
Jq is the zero order Bessel function and F{z) is a function 
that describes the redshift evolution of density fiuctua- 
tions and any evolution of th e galaxy population. Fo\- 
lowing lPodelson et alJ l|2002|) we set F(z) = 1, which is 
a much better approximation for the very local 2MASS 
data than for other surveys. ^ is the probability distri- 
bution of galaxy redshifts in the survey, i.e. the number 
of galaxies per redshift bin normalized to unity. The co- 
moving distance is 



dz' 



(13) 
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Fig. 15. — Relative redshift distributions of 2MASS galaxies 
identified in the SDSS EDR (red, N = 6238) and the 2dFGRS 100k 
release (blue, N = 9649). The total distribution for the 15,887 
sources is given in bold with Poisson errors shown for each bin. The 
smooth curve is our best fit to equation (15), which has median 
redshift Zm = 0.0741. 



where 



Eiz): 



dz 
dX 



(14) 



Van(l+2)3 + r!if(l+z)2 + r!A 



ijPeebleallQSOfl . Thus, with the measured angular corre- 
lation function and knowledge of the redshift distribution 
of the galaxies in the survey, it is possible to estimate 
P(k). We adopt the currently fashionable ACDM model 
(jSpergel et al.l l200^. ri,„ = 0.3 and JIa — 0.7, and note 
that the depth of 2MASS makes cosmological dependence 
small. 

There are a number of methods for performing the 
inversion. We will use the method of Singular Value 
Decomposition advocated by 'Eisenstei n fc Zalda rriaga' 
200 1.), and used for the SDSS analysis of the EDR 



Jd 



Dodelson 

5.1. Redshift Distribution 

To determine the redshift distribution of 2MASS 
galaxies we use redshifts measured by the SDSS EDR 
i Stoughton et a l. 2002) and by the 2dFGRS 100k release 
i Colless et aOlloOlli . We identify 6,238 2MASS galax- 
ies with redshifts in the SDSS EDR and 9, 649 2MASS 
galaxies with redshifts in the 2dFGRS 100k release. The 
distribution of redshifts from the two samples (shown 
in Figure I15|l are similar, although statistically signifi- 
cant differences, which are caused by the local large-scale 
structure, are apparent at most redshifts. 

For the combined sample of 15, 887 redshifts we fit the 
functional form 



dN _ 3-2^ 
H ~ 2{zmllAl2f 



expl 



-(1.412z/z„)3/2] (15) 
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Fig. 16. — The absolute value of the kernel as a function of kO. 
The dashed line shows where the kernel is negative. The inversion 
will be most sensitive to the first zero of the kernel, which occurs 
at ke = 0.013. 



(jBaugh fc Efstat hiou'1993D. where is the median red- 
shift. We use bootstrap resampling to estimate the 
uncertainty in our determination of and find that 
Zm = 0.0741 ±0.004. This statistical uncertainty is prob- 
ably a slight underestimate of the cosmic variance uncer- 
tainty as evidenced from the median redshifts of the two 
subsamples, which is 0.0746 and 0.0735 for 2dFGRS and 
SDSS, respectively. These are both more than one sigma 
from the median of the combined sample. We compare 
the functional form dN/dz to the data in Figure [TKl The 
kernel resulting from using this dN / dz in equation El is 
shown in Figure ITHl 

5.2. Inverting the Angular Correlation Function 

In practice we measure the angular correlation function 
in Ng bins and would like to estimate the power spectrum 
in Nk bins. We can transform the integral in equation 
Hll(l into a discrete sum and write the equation in matrix 
form as 



(16) 



where w is an Nq element vector containing the value of 
w{9) in each bin, P is an element vector containing 
the values of P{k) in each bin, and the Ng x Nk matrix 
G is a discretization of equation ((TT]) . 

The inversion is thus reduced to a standard problem in 
linear algebra; the best fit power spectrum being given by 
P — G~^w and the power spectrum covariance matrix by 
Cp = G*C^„ G- Singular Value (SV) Decomposition (for 
a review see iPress et al.lll992l §2.6 and §15.6) provides 
an effective method of solving and analyzing this prob- 
lem. The matrix G can be decomposed into G = UWV* 
where is a square diagonal matrix of the singular val- 
ues and U and V are unitary. Then G~^ is simply given 
by G^^ — VW~^U*. However, those elements of W that 
are singular or very small will make a large contribution 
to G~^. To remove this contribution, these eigenvalues 
are set to zero in the (pseudo-) inverse matrix W~^. In 
practice, we remove eigenvalues that cause large oscilla- 
tions in the resulting power spectrum. 
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Fig. 17. — The recovered power spectrum for different choices in 
the number of SV modes retained for the inversion. Retaining too 
many modes leads to wild oscillations in P{k). With 16 SV modes 
we are sensitive to small k but do not exhibit oscillations. 



Our 9 bins span the range 0.00014 — 0.3 radians. Since 
the first zero of the kernel occurs at k9 = 0.013, we ex- 
pect sensitivity to k values between 0.04 — lOO/i Mpc~^. 
However, since small k are the most interesting for cos- 
mology as they probe the linear regime, we attempt to 
extend our analysis down to fc = 0.004. We use 18 
bins in wavenumbcr, 4 per decade, spanning the range 
0.003 — lOO/i Mpc^^. Figure El shows the recovered 
power spectrum as a function of the number of modes 
we retain from the SV decomposition. We find that in- 
cluding 16 SV modes gives us sensitivity to small k but 
doesn't result in wild oscillations of the power spectrum. 
This produces an acceptable fit to w{9) with = 4.1 
with Ng — Nsv — 9 degrees of freedom. 

In Figure 1181 we plot the inferred power spectrum 
(with 16 SV modes ) compared to the p ower spectra mea- 
sured for SDSS ( Dodelson et alJl2002L r-b and selected) 
and APM (,Eisenstein fc ZaldarriaEfal 120011 6j-band se- 
lected). There is general agreement for the power spec- 
tra measured in all three surveys. The drop in power at 
k < 0.03/i Mpc~^ is not statistically significant as demon- 
strated below. The plotted error bars are the square root 
of the diagonal elements of the inverse covariance matrix, 

\J \i ■ Since the k bins are not independent, they rep- 
resent the uncertainty in a bin if the values of all other 
bins are kept fixed. They are shown only for comparison 
to the other surveys and to get a feel for the relative er- 
ror in a given bin. To treat properly the uncertainty in 
the measured Pik) one must use the full covariance ma- 
trix. In Table 13 we give the inverse correlation matrix, 
which can be used to recover the inverse covariance ma- 
trix r~^ (contact the authors for a more precise table). 
The elements of the inverse correlation matrix are given 



by 



(17) 



The correlation matrix is shown in Figure lT^ The impor- 
tance of off diagonal elements and the oscillatory nature 
of different k bins is clearly seen. 



Fig. 18. — The power spectrum measured from 2MASS (solid 
line). For comparison the results from the APM survey (triangles) 
as analyzed by Eisenstein & Zaldarriaga 12001), and for the r = 
18 — 19 magnitude bin of the SDSS survey (squares) are shown. In 
all cases error bars are the square root of the diagonal elements of 

the inverse covariance matrix, Cp\^■ The dashed line is the best 

fit CDM type power spectrum to the linear regime k < 0.15. The 
dotted line is the power spectrum from the CMB as determined by 
ISpergel et al...200a'l . 



5.3. CDM Models 

In the CDM paradigm for the evolution of the Uni- 
verse, the power spectrum can be predicted based on 
the value of the cosmological parameters of which the 
most important are the matter density of the universe 
riM and the Hubble constant lOO/i Mpc"""^ kms""'^. The 
shape of the power spectrum is most sensitive to the 
combination D,Mh. Along with a normalization, in terms 
of cTg for example, this specifies a power spectrum given 
an initial spectral index rig. Thus iluh and tJg serve 
as a convenient parameterizations of the power spec- 
trum on linear scales. We fit P{k) inverted from 'w{9) 
only for bins with k < 0.1 5/t Mpc~^, bins that are 
still in the linear regime (see iPercival et al. 2001, Fig- 
ur e 4). We use the tran sfer function fitting formula 
of lEisenstein fc Hul l|1999D . Ug = 1.0 and take the the 
baryon density = 0.02fe~^ in agreement with re- 
sults from BBN l|Kirkman et al . 2003) and the CMB 
fSpe.rgel et al. 2003) We plot the best-fit power spectrum 
in Figure 1181 including the error in the median redshift 
of the sample. The CDM power spectrum is a good fit 
to the data giving x^/d.o.f. = 0.6. Hence, the apparent 
drop in power on large sca l es, as seen by some authors 
fCaz tanaga & BaughllloM lAllgood et al.ll200l(l . is not 
statistically significant. 

The constraints on erg and SIm^ are shown in Figure 
EH We find erg = 1.0 ± 0.1 and flnh = 0.13 ± 0.07 (95% 
confidence limits). One must remember that this is only 
a parameterized fit to the power spectrum of Kg < 13.57 
galaxies; any relation to cosmology requires an under- 
standing of how galaxies trace the underlying dark mat- 
ter density field. The 2dfGRS find a power spectrum best 
fit by O^/i = 0.2 ± 0.03 (Pcrcival ct al. 2001) consistent 
with our measurement. 

In comparison, cosmological parameters measured 
from the WMAP sateUite alone give flMh = 0.19 ± 0.03, 
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Fig. 19. — The correlation matrix of the power spectrum mea- 
sured from 2MASS. The correlation matrix is the covariance matrix 
normalized by its diagonal elements. One sees how different k bins 
are not at all independent but are all interwoven in an oscillatory 
fashion. The elements of the correlation matrix are also given in 
Table El 



erg = 0.9 ± 0.1 and = .99 llSoergel et alJl2003|) . This 
would imply a i^Tg-band linear bias of 6^ = 1.1 ± 0.2 
in good agreement with the iir<,-band bias determined 
from the 2MASS c lustering dipole of Ik = 1.37 ± 0.34 
l|Maller et al.ll2(iol . 

6. CONCLUSIONS 

We have presented a measurement of the angular cor- 
relation function and three-dimensional power spectrum 
for galaxies from the 2MASS catalog. We have mini- 
mized the contribution of the possible contaminants of 
stars, dust, seeing and sky brightness by studying their 
cross-correlation with the galaxy density and making 
cuts in the data until these cross-correlations are less 
then the auto-correlation signal. These restrictions on 
the data limit us to |6| > 20° and dust extinctions 
I^Ks < 0.05. More than a half million galaxies remain 
for estimating the angular correlation function. 

The best fit power law to the measured angular corre- 
lation function has a slope of —0.79 ± 0.02 and an am- 
plitude at one degree of 0.10 ± 0.01 out to 2.5°. There 
are oscillations around this power law that are statisti- 
cally significant. The largest oscillation occurs at about 
0.8 degrees, corresponding to 7QQh~^ kpc at the median 
redshift of our survey, as expected in halo occupation dis- 
tribution descriptions of galaxy cluste ring (Seliak 2000; 
iBerlind fc Weinberg 2002VZ ehavi et al.l2004lp . The slope 
of this power law is in good agreement with other deter- 
mina tions of w{9) (Maddo x et alJlT99(Tl: iConnollv et al.l 
l2002j) . We divide the sample into three magnitude bins 
and estimate the angular correlation function for each 
magnitude bin. We confirm that differences in the w{9) 
amplitude at these three brightnesses is fully explained 
by differences in the median redshift, consistent with true 



clustering. We partition the data by northern and south- 
ern Galactic latitude and see that at large angles the 
northern Galactic latitudes show a greater correlation 
amplitude, caused by the observed large structures of 
the Virgo super cluster and the Shapley concentration 
in th e northern Galactic hemisphere (see IMaller et alJ 
120031 Fig. 2). We also partition the data by effective sur- 
face brightness and find that galaxies with higher surface 
brightness are more clustered, which is a manifestation of 
the morphology-density relation (Dressier 1980). Their 
w{9) may also have a slightly steeper slope and might be 
evidence for nonlocal biasing. 

We solve for the three-dimensional power spectrum 
from the angular correlation function using Singular 
Value Decomposition. Our resulting power spectrum 
is in good agreement wit h other m easurements us- 
ing the same metho d (Eisenstein fc Zaldarriaga 200^ 
iDodelson et al.ll2002() . The best fit CDM power spec- 
tra gives values oi^luh = 0.13±0.07 and trg — l.OiO.l, 
assuming a spectral index of 1.0. This fit is to the galaxy 
power spectrum and its relation to cosmological param- 
eters will depend on how galaxies trace the underlying 
mass distribution. The ra tio of our measured a% to that 
determined from WMAP ijSpergel et al.ll200^ would im- 
ply a Kg band linear bias oihK = 1.1 ± 0.2. 

Our i^s-band selected power spectrum, is in good 
agreement with other determinations of the power spec- 
trum in surveys selected in other wavebands in the lin- 
ear regime. On smaller length scales there is a sig- 
nificant difference between other galaxy power spectra 
and that measured for i^s-band selected galaxies. Thus 
the power spectrum measured here, in combination with 
the power spectrum measured from the CMB and other 
large galaxy surveys, will enable one to place strong con- 
straints on theories of galaxy formation. 
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TABLE 3 

The results of the best fit power spectrum using 16 SV modes. The first row are the centers of the k bins in h Mpc"'^. The second row are the best fit values of 
P(fe). The third row is the square root of the covariance matrix, ^fCpJi- The fourth row is the square root of the inverse covariance matrix C^\^ and the rest 

OF THE TABLE GIVES THE 18 X 18 SYMMETRIC INVERSE CORRELATION MATRIX (SEE TEXT). NOTE REDSHIFT ERRORS ARE NOT INCLUDED IN THE ABOVE. FOR A TABLE WITH GREATER 

precision PLEASE CONTACT THE AUTHORS. 



